Structure, function, and phylogeny of the mating locus in the Rhizopus oryzae complex.

The Rhizopus oryzae species complex is a group of zygomycete fungi that are common, cosmopolitan saprotrophs. Some strains are used beneficially for production of Asian fermented foods but they can also act as opportunistic human pathogens. Although R. oryzae reportedly has a heterothallic (+/-) mating system, most strains have not been observed to undergo sexual reproduction and the genetic structure of its mating locus has not been characterized. Here we report on the mating behavior and genetic structure of the mating locus for 54 isolates of the R. oryzae complex. All 54 strains have a mating locus similar in overall organization to Phycomyces blakesleeanus and Mucor circinelloides (Mucoromycotina, Zygomycota). In all of these fungi, the minus (-) allele features the SexM high mobility group (HMG) gene flanked by an RNA helicase gene and a TP transporter gene (TPT). Within the R. oryzae complex, the plus (+) mating allele includes an inserted region that codes for a BTB/POZ domain gene and the SexP HMG gene. Phylogenetic analyses of multiple genes, including the mating loci (HMG, TPT, RNA helicase), ITS1-5.8S-ITS2 rDNA, RPB2, and LDH genes, identified two distinct groups of strains. These correspond to previously described sibling species R. oryzae sensu stricto and R. delemar. Within each species, discordant gene phylogenies among multiple loci suggest an outcrossing population structure. The hypothesis of random-mating is also supported by a 50:50 ratio of plus and minus mating types in both cryptic species. When crossed with tester strains of the opposite mating type, most isolates of R. delemar failed to produce zygospores, while isolates of R. oryzae produced sterile zygospores. In spite of the reluctance of most strains to mate in vitro, the conserved sex locus structure and evidence for outcrossing suggest that a normal sexual cycle occurs in both species.

Morphological and genomic characterization of Filobasidiella depauperata: a homothallic sibling species of the pathogenic cryptococcus species complex.

The fungal species Cryptococcus neoformans and Cryptococcus gattii cause respiratory and neurological disease in animals and humans following inhalation of basidiospores or desiccated yeast cells from the environment. Sexual reproduction in C. neoformans and C. gattii is controlled by a bipolar system in which a single mating type locus (MAT) specifies compatibility. These two species are dimorphic, growing as yeast in the asexual stage, and producing hyphae, basidia, and basidiospores during the sexual stage. In contrast, Filobasidiella depauperata, one of the closest related species, grows exclusively as hyphae and it is found in association with decaying insects. Examination of two available strains of F. depauperata showed that the life cycle of this fungal species shares features associated with the unisexual or same-sex mating cycle in C. neoformans. Therefore, F. depauperata may represent a homothallic and possibly an obligately sexual fungal species. RAPD genotyping of 39 randomly isolated progeny from isolate CBS7855 revealed a new genotype pattern in one of the isolated basidiospores progeny, therefore suggesting that the homothallic cycle in F. depauperata could lead to the emergence of new genotypes. Phylogenetic analyses of genes linked to MAT in C. neoformans indicated that two of these genes in F. depauperata, MYO2 and STE20, appear to form a monophyletic clade with the MATa alleles of C. neoformans and C. gattii, and thus these genes may have been recruited to the MAT locus before F. depauperata diverged. Furthermore, the ancestral MATa locus may have undergone accelerated evolution prior to the divergence of the pathogenic Cryptococcus species since several of the genes linked to the MATa locus appear to have a higher number of changes and substitutions than their MATalpha counterparts. Synteny analyses between C. neoformans and F. depauperata showed that genomic regions on other chromosomes displayed conserved gene order. In contrast, the genes linked to the MAT locus of C. neoformans showed a higher number of chromosomal translocations in the genome of F. depauperata. We therefore propose that chromosomal rearrangements appear to be a major force driving speciation and sexual divergence in these closely related pathogenic and saprobic species.

Evolution of the sex-related locus and genomic features shared in microsporidia and fungi.

BACKGROUND: Microsporidia are obligate intracellular, eukaryotic pathogens that infect a wide range of animals from nematodes to humans, and in some cases, protists. The preponderance of evidence as to the origin of the microsporidia reveals a close relationship with the fungi, either within the kingdom or as a sister group to it. Recent phylogenetic studies and gene order analysis suggest that microsporidia share a particularly close evolutionary relationship with the zygomycetes. METHODOLOGY/PRINCIPAL FINDINGS: Here we expanded this analysis and also examined a putative sex-locus for variability between microsporidian populations. Whole genome inspection reveals a unique syntenic gene pair (RPS9-RPL21) present in the vast majority of fungi and the microsporidians but not in other eukaryotic lineages. Two other unique gene fusions (glutamyl-prolyl tRNA synthetase and ubiquitin-ribosomal subunit S30) that are present in metazoans, choanoflagellates, and filasterean opisthokonts are unfused in the fungi and microsporidians. One locus previously found to be conserved in many microsporidian genomes is similar to the sex locus of zygomycetes in gene order and architecture. Both sex-related and sex loci harbor TPT, HMG, and RNA helicase genes forming a syntenic gene cluster. We sequenced and analyzed the sex-related locus in 11 different Encephalitozoon cuniculi isolates and the sibling species E. intestinalis (3 isolates) and E. hellem (1 isolate). There was no evidence for an idiomorphic sex-related locus in this Encephalitozoon species sample. According to sequence-based phylogenetic analyses, the TPT and RNA helicase genes flanking the HMG genes are paralogous rather than orthologous between zygomycetes and microsporidians. CONCLUSION/SIGNIFICANCE: The unique genomic hallmarks between microsporidia and fungi are independent of sequence based phylogenetic comparisons and further contribute to define the borders of the fungal kingdom and support the classification of microsporidia as unusual derived fungi. And the sex/sex-related loci appear to have been subject to frequent gene conversion and translocations in microsporidia and zygomycetes.

Comparative analyses of clinical and environmental populations of Cryptococcus neoformans in Botswana.

Cryptococcus neoformans var. grubii (Cng) is the most common cause of fungal meningitis, and its prevalence is highest in sub-Saharan Africa. Patients become infected by inhaling airborne spores or desiccated yeast cells from the environment, where the fungus thrives in avian droppings, trees and soil. To investigate the prevalence and population structure of Cng in southern Africa, we analysed isolates from 77 environmental samples and 64 patients. We detected significant genetic diversity among isolates and strong evidence of geographic structure at the local level. High proportions of isolates with the rare MATa allele were observed in both clinical and environmental isolates; however, the mating-type alleles were unevenly distributed among different subpopulations. Nearly equal proportions of the MATa and MATα mating types were observed among all clinical isolates and in one environmental subpopulation from the eastern part of Botswana. As previously reported, there was evidence of both clonality and recombination in different geographic areas. These results provide a foundation for subsequent genomewide association studies to identify genes and genotypes linked to pathogenicity in humans.

Environmental and genetic determinants of colony morphology in yeast.

Nutrient stresses trigger a variety of developmental switches in the budding yeast Saccharomyces cerevisiae. One of the least understood of such responses is the development of complex colony morphology, characterized by intricate, organized, and strain-specific patterns of colony growth and architecture. The genetic bases of this phenotype and the key environmental signals involved in its induction have heretofore remained poorly understood. By surveying multiple strain backgrounds and a large number of growth conditions, we show that limitation for fermentable carbon sources coupled with a rich nitrogen source is the primary trigger for the colony morphology response in budding yeast. Using knockout mutants and transposon-mediated mutagenesis, we demonstrate that two key signaling networks regulating this response are the filamentous growth MAP kinase cascade and the Ras-cAMP-PKA pathway. We further show synergistic epistasis between Rim15, a kinase involved in integration of nutrient signals, and other genes in these pathways. Ploidy, mating-type, and genotype-by-environment interactions also appear to play a role in the controlling colony morphology. Our study highlights the high degree of network reuse in this model eukaryote; yeast use the same core signaling pathways in multiple contexts to integrate information about environmental and physiological states and generate diverse developmental outputs.

Cryptococcal cell morphology affects host cell interactions and pathogenicity.

Cryptococcus neoformans is a common life-threatening human fungal pathogen. The size of cryptococcal cells is typically 5 to 10 microm. Cell enlargement was observed in vivo, producing cells up to 100 microm. These morphological changes in cell size affected pathogenicity via reducing phagocytosis by host mononuclear cells, increasing resistance to oxidative and nitrosative stress, and correlated with reduced penetration of the central nervous system. Cell enlargement was stimulated by coinfection with strains of opposite mating type, and ste3aDelta pheromone receptor mutant strains had reduced cell enlargement. Finally, analysis of DNA content in this novel cell type revealed that these enlarged cells were polyploid, uninucleate, and produced daughter cells in vivo. These results describe a novel mechanism by which C. neoformans evades host phagocytosis to allow survival of a subset of the population at early stages of infection. Thus, morphological changes play unique and specialized roles during infection.

Similar T-cell immune responses induced by group M consensus env immunogens with wild-type or minimum consensus variable regions.

Consensus HIV-1 genes can decrease the genetic distances between candidate immunogens and field virus strains. To ensure the functionality and optimal presentation of immunologic epitopes, we generated two group-M consensus env genes that contain variable regions either from a wild-type B/C recombinant virus isolate (CON6) or minimal consensus elements (CON-S) in the V1, V2, V4, and V5 regions. C57BL/6 and BALB/c mice were primed twice with CON6, CON-S, and subtype control (92UG37_A and HXB2/Bal_B) DNA and boosted with recombinant vaccinia virus (rVV). Mean antibody titers against 92UG37_A, 89.6_B, 96ZM651_C, CON6, and CON-S Env protein were determined. Both CON6 and CON-S induced higher mean antibody titers against several of the proteins, as compared with the subtype controls. However, no significant differences were found in mean antibody titers in animals immunized with CON6 or CON-S. Cellular immune responses were measured by using five complete Env overlapping peptide sets: subtype A (92UG37_A), subtype B (MN_B, 89.6_B and SF162_B), and subtype C (Chn19_C). The intensity of the induced cellular responses was measured by using pooled Env peptides; T-cell epitopes were identified by using matrix peptide pools and individual peptides. No significant differences in T-cell immune-response intensities were noted between CON6 and CON-S immunized BALB/c and C57BL/6 mice. In BALB/c mice, 10 and eight nonoverlapping T-cell epitopes were identified in CON6 and CON-S, whereas eight epitopes were identified in 92UG37_A and HXB2/BAL_B. In C57BL/6 mice, nine and six nonoverlapping T-cell epitopes were identified after immunization with CON6 and CON-S, respectively, whereas only four and three were identified in 92UG37_A and HXB2/BAL_B, respectively. When combined together from both mouse strains, 18 epitopes were identified. The group M artificial consensus env genes, CON6 and CON-S, were equally immunogenic in breadth and intensity for inducing humoral and cellular immune responses.

Recursive directional ligation by plasmid reconstruction allows rapid and seamless cloning of oligomeric genes.

This paper reports a new strategy, recursive directional ligation by plasmid reconstruction (PRe-RDL), to rapidly clone highly repetitive polypeptides of any sequence and specified length over a large range of molecular weights. In a single cycle of PRe-RDL, two halves of a parent plasmid, each containing a copy of an oligomer, are ligated together, thereby dimerizing the oligomer and reconstituting a functional plasmid. This process is carried out recursively to assemble an oligomeric gene with the desired number of repeats. PRe-RDL has several unique features that stem from the use of type IIs restriction endonucleases: first, PRe-RDL is a seamless cloning method that leaves no extraneous nucleotides at the ligation junction. Because it uses type IIs endonucleases to ligate the two halves of the plasmid, PRe-RDL also addresses the major limitation of RDL in that it abolishes any restriction on the gene sequence that can be oligomerized. The reconstitution of a functional plasmid only upon successful ligation in PRe-RDL also addresses two other limitations of RDL: the significant background from self-ligation of the vector observed in RDL, and the decreased efficiency of ligation due to nonproductive circularization of the insert. PRe-RDL can also be used to assemble genes that encode different sequences in a predetermined order to encode block copolymers or append leader and trailer peptide sequences to the oligomerized gene.

Chiropteran types I and II interferon genes inferred from genome sequencing traces by a statistical gene-family assembler.

BACKGROUND: The rate of emergence of human pathogens is steadily increasing; most of these novel agents originate in wildlife. Bats, remarkably, are the natural reservoirs of many of the most pathogenic viruses in humans. There are two bat genome projects currently underway, a circumstance that promises to speed the discovery host factors important in the coevolution of bats with their viruses. These genomes, however, are not yet assembled and one of them will provide only low coverage, making the inference of most genes of immunological interest error-prone. Many more wildlife genome projects are underway and intend to provide only shallow coverage. RESULTS: We have developed a statistical method for the assembly of gene families from partial genomes. The method takes full advantage of the quality scores generated by base-calling software, incorporating them into a complete probabilistic error model, to overcome the limitation inherent in the inference of gene family members from partial sequence information. We validated the method by inferring the human IFNA genes from the genome trace archives, and used it to infer 61 type-I interferon genes, and single type-II interferon genes in the bats Pteropus vampyrus and Myotis lucifugus. We confirmed our inferences by direct cloning and sequencing of IFNA, IFNB, IFND, and IFNK in P. vampyrus, and by demonstrating transcription of some of the inferred genes by known interferon-inducing stimuli. CONCLUSION: The statistical trace assembler described here provides a reliable method for extracting information from the many available and forthcoming partial or shallow genome sequencing projects, thereby facilitating the study of a wider variety of organisms with ecological and biomedical significance to humans than would otherwise be possible.

A genome-wide association study identifies susceptibility variants for type 2 diabetes in Han Chinese.

To investigate the underlying mechanisms of T2D pathogenesis, we looked for diabetes susceptibility genes that increase the risk of type 2 diabetes (T2D) in a Han Chinese population. A two-stage genome-wide association (GWA) study was conducted, in which 995 patients and 894 controls were genotyped using the Illumina HumanHap550-Duo BeadChip for the first genome scan stage. This was further replicated in 1,803 patients and 1,473 controls in stage 2. We found two loci not previously associated with diabetes susceptibility in and around the genes protein tyrosine phosphatase receptor type D (PTPRD) (P = 8.54x10(-10); odds ratio [OR] = 1.57; 95% confidence interval [CI] = 1.36-1.82), and serine racemase (SRR) (P = 3.06x10(-9); OR = 1.28; 95% CI = 1.18-1.39). We also confirmed that variants in KCNQ1 were associated with T2D risk, with the strongest signal at rs2237895 (P = 9.65x10(-10); OR = 1.29, 95% CI = 1.19-1.40). By identifying two novel genetic susceptibility loci in a Han Chinese population and confirming the involvement of KCNQ1, which was previously reported to be associated with T2D in Japanese and European descent populations, our results may lead to a better understanding of differences in the molecular pathogenesis of T2D among various populations.

G-protein-coupled receptor genes as protooncogenes: constitutively activating mutation of the alpha 1B-adrenergic receptor enhances mitogenesis and tumorigenicity.

The alpha 1B-adrenergic receptor (alpha 1B-ADR) is a member of the G-protein-coupled family of transmembrane receptors. When transfected into Rat-1 and NIH 3T3 fibroblasts, this receptor induces focus formation in an agonist-dependent manner. Focus-derived, transformed fibroblasts exhibit high levels of functional alpha 1B-ADR expression, demonstrate a catecholamine-induced enhancement in the rate of cellular proliferation, and are tumorigenic when injected into nude mice. Induction of neoplastic transformation by the alpha 1B-ADR, therefore, identifies this normal cellular gene as a protooncogene. Mutational alteration of this receptor can lead to activation of this protooncogene, resulting in an enhanced ability of agonist to induce focus formation with a decreased latency and quantitative increase in transformed foci. In contrast to cells expressing the wild-type alpha 1B-ADR, focus formation in "oncomutant"-expressing cell lines appears constitutively activated with the generation of foci in unstimulated cells. Further, these cell lines exhibit near-maximal rates of proliferation even in the absence of catecholamine supplementation. They also demonstrate an enhanced ability for tumor generation in nude mice with a decreased period of latency compared with cells expressing the wild-type receptor. Thus, the alpha 1B-ADR gene can, when overexpressed and activated, function as an oncogene inducing neoplastic transformation. Mutational alteration of this receptor gene can result in the activation of this protooncogene, enhancing its oncogenic potential. These findings suggest that analogous spontaneously occurring mutations in this class of receptor proteins could play a key role in the induction or progression of neoplastic transformation and atherosclerosis.

Analysis of the genome and transcriptome of Cryptococcus neoformans var. grubii reveals complex RNA expression and microevolution leading to virulence attenuation.

Cryptococcus neoformans is a pathogenic basidiomycetous yeast responsible for more than 600,000 deaths each year. It occurs as two serotypes (A and D) representing two varieties (i.e. grubii and neoformans, respectively). Here, we sequenced the genome and performed an RNA-Seq-based analysis of the C. neoformans var. grubii transcriptome structure. We determined the chromosomal locations, analyzed the sequence/structural features of the centromeres, and identified origins of replication. The genome was annotated based on automated and manual curation. More than 40,000 introns populating more than 99% of the expressed genes were identified. Although most of these introns are located in the coding DNA sequences (CDS), over 2,000 introns in the untranslated regions (UTRs) were also identified. Poly(A)-containing reads were employed to locate the polyadenylation sites of more than 80% of the genes. Examination of the sequences around these sites revealed a new poly(A)-site-associated motif (AUGHAH). In addition, 1,197 miscRNAs were identified. These miscRNAs can be spliced and/or polyadenylated, but do not appear to have obvious coding capacities. Finally, this genome sequence enabled a comparative analysis of strain H99 variants obtained after laboratory passage. The spectrum of mutations identified provides insights into the genetics underlying the micro-evolution of a laboratory strain, and identifies mutations involved in stress responses, mating efficiency, and virulence.

Different Mechanisms Confer Gradual Control and Memory at Nutrient- and Stress-Regulated Genes in Yeast.

Cells respond to environmental stimuli by fine-tuned regulation of gene expression. Here we investigated the dose-dependent modulation of gene expression at high temporal resolution in response to nutrient and stress signals in yeast. The GAL1 activity in cell populations is modulated in a well-defined range of galactose concentrations, correlating with a dynamic change of histone remodeling and RNA polymerase II (RNAPII) association. This behavior is the result of a heterogeneous induction delay caused by decreasing inducer concentrations across the population. Chromatin remodeling appears to be the basis for the dynamic GAL1 expression, because mutants with impaired histone dynamics show severely truncated dose-response profiles. In contrast, the GRE2 promoter operates like a rapid off/on switch in response to increasing osmotic stress, with almost constant expression rates and exclusively temporal regulation of histone remodeling and RNAPII occupancy. The Gal3 inducer and the Hog1 mitogen-activated protein (MAP) kinase seem to determine the different dose-response strategies at the two promoters. Accordingly, GAL1 becomes highly sensitive and dose independent if previously stimulated because of residual Gal3 levels, whereas GRE2 expression diminishes upon repeated stimulation due to acquired stress resistance. Our analysis reveals important differences in the way dynamic signals create dose-sensitive gene expression outputs.

Characterization of the FKBP12-Encoding Genes in Aspergillus fumigatus.

Invasive aspergillosis, largely caused by Aspergillus fumigatus, is responsible for a growing number of deaths among immunosuppressed patients. Immunosuppressants such as FK506 (tacrolimus) that target calcineurin have shown promise for antifungal drug development. FK506-binding proteins (FKBPs) form a complex with calcineurin in the presence of FK506 (FKBP12-FK506) and inhibit calcineurin activity. Research on FKBPs in fungi is limited, and none of the FKBPs have been previously characterized in A. fumigatus. We identified four orthologous genes of FKBP12, the human FK506 binding partner, in A. fumigatus and designated them fkbp12-1, fkbp12-2, fkbp12-3, and fkbp12-4. Deletional analysis of the four genes revealed that the Δfkbp12-1 strain was resistant to FK506, indicating FKBP12-1 as the key mediator of FK506-binding to calcineurin. The endogenously expressed FKBP12-1-EGFP fusion protein localized to the cytoplasm and nuclei under normal growth conditions but also to the hyphal septa following FK506 treatment, revealing its interaction with calcineurin. The FKBP12-1-EGFP fusion protein didn't localize at the septa in the presence of FK506 in the cnaA deletion background, confirming its interaction with calcineurin. Testing of all deletion strains in the Galleria mellonella model of aspergillosis suggested that these proteins don't play an important role in virulence. While the Δfkbp12-2 and Δfkbp12-3 strains didn't show any discernable phenotype, the Δfkbp12-4 strain displayed slight growth defect under normal growth conditions and inhibition of the caspofungin-mediated "paradoxical growth effect" at higher concentrations of the antifungal caspofungin. Together, these results indicate that while only FKBP12-1 is the bona fide binding partner of FK506, leading to the inhibition of calcineurin in A. fumigatus, FKBP12-4 may play a role in basal growth and the caspofungin-mediated paradoxical growth response. Exploitation of differences between A. fumigatus FKBP12-1 and human FKBP12 will be critical for the generation of fungal-specific FK506 analogs to inhibit fungal calcineurin and treat invasive fungal disease.

Pathogen-Specific Adaptations to Conserved Signaling Pathways in Cryptococcus neoformans

Cryptococcus neoformans is an opportunistic fungal pathogen that causes significant disease worldwide. Even though this fungus has not evolved specifically to cause human disease, it has a remarkable ability to adapt to many different environments within its infected host. C. neoformans adapts by utilizing conserved eukaryotic and fungal-specific signaling pathways to sense and respond to stresses within the host. Upon infection, two of the most significant environmental changes this organism experiences are elevated temperature and high pH.

Conserved Rho and Ras family GTPases are central regulators of thermotolerance in C. neoformans. Many GTPases require prenylation to associate with cellular membranes and function properly. Using molecular genetic techniques, microscopy, and infection models, I demonstrated that the prenyltransferase, geranylgeranyl transferase I (GGTase I) is required for thermotolerance and pathogenesis. Using fluorescence microscopy, I found that only a subset of conserved GGTase I substrates requires this enzyme for membrane localization. Therefore, the C. neoformans GGTase I may recognize its substrate in a slightly different manner than other eukaryotic organisms.

The alkaline response transcription factor, Rim101, is a central regulator of stress-response genes important for adapting to the host environment. In particular, Rim101 regulates cell surface alterations involved in immune avoidance. In other fungi, Rim101 is activated by alkaline pH through a conserved signaling pathway, but this pathway had yet been characterized in C. neoformans. Using molecular genetic techniques, I identified and analyzed the conserved members of the Rim pathway. I found that it was only partially conserved in C. neoformans, missing the components that sense pH and initiate pathway activation. Using a genetic screen, I identified a novel Rim pathway component named Rra1. Structural prediction and genetic epistasis experiments suggest that Rra1 may serve as the Rim pathway pH sensor in C. neoformans and other related basidiomycete fungi.

To explore the relevance of Rim pathway signaling in the interaction of C neoformans with its host, I characterized the Rim101-regulated cell wall changes that prevent immune detection. Using HPLC, enzymatic degradation, and cell wall stains, I found that the rim101Δ mutation resulted in increased cell wall chitin exposure. In vitro co-culture assays demonstrated that increased chitin exposure is associated with enhanced activation of macrophages and dendritic cells. To further test this association, I demonstrated that other mutant strains with increased chitin exposure induce macrophage and dendritic cell responses similar to rim101Δ. We used primary macrophages from mutant mouse lines to demonstrate that members of both the Toll-like receptor and C-type lectin receptor families are involved in detecting strains with increased chitin exposure. Finally, in vivo immunological experiments demonstrated that the rim101Δ strain induced a global inflammatory immune response in infected mouse lungs, expanding upon our previous in vivo rim101Δ studies. These results demonstrate that cell wall organization largely determines how fungal cells are detected by the immune system.